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Abstract. We review recent results, obtained with P. Fendley, on frustration of quantum charges in lattice 
models for itinerant fermions with strong repulsive interactions. A judicious tuning of kinetic and interac- 
tion terms leads to models possessing supersymmetry. In such models frustration takes the form of what 
we call superfrustration; an extensive degeneracy of supersymmetric ground states. We present a gallery 
of examples of superfrustration on a variety of 2D lattices. 
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1 Introduction 

When charged particles, with repulsive interactions, are 
placed on a lattice one expects geometric frustration: 
depending on the lattice and the number of particles, there 
can be many configurations that realize the lowest possi- 
ble interaction energy. Well-studied examples are charges 
on the triangular and checkerboard lattices (see [T] for 
some early references). In general, including kinetic terms 
for the quantum charges lifts the degeneracies. Depending 
on details, this may give rise to novel phases of quantum 
matter with remarkable physical properties |2j- The the- 
oretical tools for studying these systems are limited: one 
typically relies on a strong coupling expansions and on 
numerics. 

Recent work by P. Fendley and one of the authors |3j 
has uncovered models for strongly interacting itinerant 
fermions which display a strong form of quantum charge 
frustration, which we call superfrustration. These mod- 
els (defined on 2D or 3D lattices) have a large, exact 
ground state degeneracy in the presence of kinetic terms. 
Superfrustration thus arises due to a subtle interplay be- 
tween kinetic terms and strong repulsive interactions. 

The term 'superfrustration' has its origin in a key prop- 
erty used to identify the models and to study their prop- 
erties, which is supersymmetry. Quite remarkably, the 
notion of supersymmetry, which was developed in the con- 
text of high energy physics, turns out to be a powerful tool 
in the analysis of strongly correlated itinerant fermions. 
This was first recognized in the context of ID models [H 
[S]. The extension to 2D and 3D then led to the discovery 
of the phenomenon of superfrustration. 

In this review, we shall first explain (Sect. [2]) how su- 
persymmetry is put to work in lattice models of correlated 



fermions. We introduce some basic tools, such as the Wit- 
ten index, make a connection to cohomology theory and 
discuss a model on a ID chain. We then move to models on 
2D lattices (Sect. [3]), where we present a heuristic geomet- 
ric intuition (3-rule) and discuss the methods employed 
in the analysis. In Sect. |3]we present a gallery of exam- 
ples, each chosen such as to illuminate specific aspects and 
features. They firmly establish the notion of superfrustra- 
tion in its various guises. We close (Sect. O with some 
thoughts about the nature of the various ground states 
and quantum phases, in particular in relation to quantum 
criticality. 



2 Supersymmetry 

2.1 Basic algebra and Hamiltonian 

In quantum mechanics, supersymmetric theories are char- 
acterized by a positive definite energy spectrum and a 
twofold degeneracy of each non-zero energy level. The two 
states with the same energy are called superpartners and 
are related by the nilpotent supercharge operator. Let us 
consider an TV = 2 supersymmetric theory, defined by two 
nilpotent supercharges Q and 
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and the Hamiltonian given by 

H = {Q\Q}. 

From this definition it follows directly that H is positive 
definite: 

= IOI^)P + |Qt|v-)|2 >o . 
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Furthermore, both Q and commute with the Hamil- 
tonian, which gives rise to the twofold degeneracy in the 
energy spectrum. In other words, all eigenstates with an 
energy Eg > form doublet representations of the su- 
persymmetry algebra. A doublet consists of two states 
\s),Q\s), such that Q^js) = 0. Finally, all states with zero 
energy must be singlets: Q\g) = Q^g) = and conversely, 
all singlets must be zero energy states [6^. In addition to 
supersymmetry our models also have a fermion-number 
symmetry generated by the operator F with 

[F, Qt] = _Qt and [F, Q] = Q. 

Consequently, F commutes with the Hamiltonian. 

The supersymmetric theories that we discuss in this 
paper describe fermionic particles. One might expect that, 
in order to be supersymmetric, these theories would need 
bosonic particles as well but this is not the case. The 
crux is that the quantum states in these theories come 
in two types: bosonic states, having an even number / of 
fermionic particles, and fermionic states with / odd. From 
the commutators of F with the supercharges, one finds 
that Q and change / by plus or minus one one unit, 
so that the supercharges map bosonic states to fermionic 
states and vice versa. 

We now make things concrete and define supersym- 
metric models for spin-less fermions on a lattice or graph 
with L sites in any dimension, following [3]). The opera- 
tor that creates a fermion on site i is written as with 
{cl,Cj} — 5ij. A simple choice for the first supercharge 
would be Q = J2i^h where the sum is over all lattice 
sites. This leads to a trivial Hamiltonian: H = L, where 
L is the number of lattice sites. To obtain a non-trivial 
Hamiltonian, we dress the fermion with a projection op- 
erator: P<j> = nj„exttoi(l ~ c]cj), which requires aU 
sites adjacent to site i to be empty. With Q = clP<i> 
and ~ X]cjP<i>, the Hamiltonian of these hard-core 
fermions reads 

H = {Q^ Q} = Y. E p<^>4c3P<3> + E p<^>- 

i j next to i i 

The first term is just a nearest neighbor hopping term 
for hard-core fermions, the second term contains a next- 
nearest neighbor repulsion, a chemical potential and a con- 
stant. The details of the latter terms will depend on the 
lattice we choose. 

Note that all the parameters in the Hamiltonian (the 
hopping t, the nearest neighbor repulsion Vi, the next- 
nearest neighbor repulsion V2 and the chemical potential 
n) are fixed by the choice of the supercharges and the 
requirement of supersymmetry and eventually the lattice. 

2.2 Witten index 

We have already discussed how supersymmetry gives rise 
to certain properties of the spectrum of the system, such 
as positivity of the energies and pairing of the excited 



states. An important issue is whether or not supersym- 
metric ground states at zero energy occur. For this one 
considers the Witten index 

W^ = tr[(-l)^e-'^^] . (1) 

Remember that all excited states come in doublets with 
the same energy and differing in their fermion-number by 
one. This means that in the trace all contributions of ex- 
cited states will cancel pairwise, and that the only states 
contributing are the zero energy ground states. We can 
thus evaluate W in the limit of /? ^ 0, where all states 
contribute (—1)^. It also follows that \ W\ is a lower bound 
to the number of zero energy ground states. 

2.3 Example: 6-site chain 

Let us consider as an example of all the above, the chain 
of six sites with periodic boundary conditions. The first 
thing we note is that the Hamiltonian for an L-site chain 
with periodic boundary conditions can be rewritten in the 
following form: 

H = i/kin + Hpot, (2) 

where 

L 

-ffkin = ^ |^Pi_i(c|ci+i -I- c\j^iCi)Pi+2 , 
i=l 
L 

Hpot — Pi~lPi+l 
i=l 
L 

= '^[1 -2ni + nini+2] 

L 

= ^{n,n,+2)+L-2F. 

i=l 

Here Pi = 1 — rii, rii = c[ci is the usual number operator 
and F = rii is the total number of fermions. (We shall 
denote eigenvalues of this operator by / and write the 
fermion density or filling fraction as = f/L.) The form 
of the hamiltonian makes clear that the hopping param- 
eter t is tuned to be equal to the next-nearest neighbor 
repulsion V2, which is tuned to unity. The nearest neigh- 
bor repulsion Vi is by definition infinite and the chemical 
potential /i is 2. Finally, there is a constant contribution 
L to the Hamiltonian. Note that the second term in the 
Hamiltonian i/pot suggests that the energy is minimized 
when the hard-core fermions are three sites apart. 

Let us consider the possible configurations of the 6- 
site chain. In addition to the empty state, there are six 
configurations with one fermion, nine with two fermions 
and two with three fermions (see Fig. [1]). Because of the 
hard-core character of the fermions, half-filling is the max- 
imal density. Clearly, the operator Q gives zero on these 
maximally filled states. On the other hand, Qt acts non- 
trivially on these states, so two of the nine states with two 
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Fig. 1. Configurations and spectrum of the 6-site chain. 

fermions are superpartners of the maximally filled states. 
The empty state |0) has an energy E = 6 and Q^O) = 0, 
whereas Q|0) = J2i4\0), so (|0),Q|0)) make up a dou- 
blet. The other five states with one fermion are annihi- 
lated by a-iid Q acts non-trivially on them, so they form 
supersymmetry doublets with five two-fermion-states. At 
this point, seven of the nine two-fermion-states are paired 
up in doublets, either with one- or three-fermion-states. 
The remaining two states cannot be part of a doublet, 
which implies that they must be singlet states and thus 
have zero energy. So we find that the 6-site chain has 
a twofold degenerate zero energy ground state at filling 
u = f/L = 1/3. The full spectrum of the 6-sitc chain is 
shown in Fig. [T] 

We observe that the ground state filling fraction of 1 /3 
agrees with the expectation that fermions tend to be three 
sites apart. This geometric rule suggests three possible 
ground states; in the full quantum theory two are realized 
as zero-energy states. Note that the actual ground state 
wavefunctions are superpositions of many different con- 
figurations. With a bit more work one can show that the 
ground states have eigenvalues exp(±7ri/3) under transla- 
tion by one site. 

Now let us compute the Witten index for this exam- 
ple. Remember that for a supersymmetric theory it simply 
reads 

W = tr(-l)^. 

Note that we can take any basis of states we like to com- 
pute the trace. Above we have specified a basis by consid- 
ering all the possible configurations of up to three fermions 
on the chain. It immediately gives M^ = l — 6 + 9 — 2 = 2 
in agreement with the existence of the two ground states 
that we found. 

We close this section with two comments. First, we 
stress that the extremely simple computation of W alone 
guarantees the existence of at least two ground states at 
zero energy. Similar results are easily established for much 
larger systems, where a direct evaluation of the ground 
state energies is way out of reach, showing the power of 
supersymmetry. Second, we observe that here the Witten 
index is exactly equal to the number of ground states. 
We will encounter examples where ground states exist at 
more than one fermion number /, leading to cancellations 



in the Witten index so that \W\ is strictly smaller than 
the number of ground states. 

2.4 Cohomology 

Supersymmetry supplies us with another tool, besides the 
Witten index, to study the ground states of the fermion 
models. This so-called cohomology method is more in- 
volved but it reveals more information about the ground 
state structure, in that it specifies the number of ground 
states for given fermion number /. 

The key ingredient is the fact that ground states are 
singlets, they are annihilated both by Q and Q^. This 
means that a ground state \g) is in the kernel of Q: Q\g) — 
0. Such a ground state is not in the image of Q, because 
if we could write \g) = Q\f), then (|/), \g)), would be a 
doublet. Equivalently, we can say that a ground state is 
closed but not exact. So the ground states span a sub- 
space Hq of the Hilbert space TC of states, such that 
Hq = kerQ/Im Q. This is precisely the definition of the 
cohomology of Q. So the ground states of a supersym- 
metric theory are in one-to-one correspondence with the 
cohomology of Q. Two states |si) and \s2) are said to be in 
the same cohomology-class if |si) — \s2) + Qls^) for some 
state I S3). Since a ground state is annihilated by both Q 
and Q^, different (i.e. linearly independent) ground states 
must be in different cohomology-classes of Q. Finally, the 
number of independent ground states is precisely the di- 
mension of the cohomology of Q and the fermion-number 
of a ground state is the same as that of the corresponding 
cohomology-class. 

There are several techniques to compute the cohomol- 
ogy, which we shall illustrate by working out examples in 
the following sections. 

2.5 Example: ID chains 

In previous work [3l|4l[5l[7| the supersymmetric model on 
the chain was studied extensively. We will summarize some 
of the results, but mostly use this case to illustrate the 
power of the tools we have developed in the previous sec- 
tions. Let us first compute the Witten index. In the exam- 
ple of the 6-site chain we saw that the Witten index can 
be computed by simply summing over all possible config- 
urations with the appropriate sign. However, because of 
the hard-core character of the fermions this is not a trivial 
problem for larger sizes. Here we shall exploit a much more 
elegant method, which will turn out very useful when we 
extend our model to more complex lattices. This method 
consists of the following steps: first divide the lattice into 
two sublattices Si and ^2. Then fix the configuration on 
Si and sum (—1)^ for the configurations on 5*2. Finally, 
sum the results over the configurations of 5*1. Of course 
the trick is to make a smart choice for the sublattices. 
For the periodic chain with L = 3j sites, we take S2 to 
be every third site. All the sites on S'2 are disconnected 
and thus every site can be either empty or occupied given 
that its neighboring sites on Si are empty. This means 



4 



Liza Huijse, Kareljan Schoutens: Superfrustration of charge degrees of freedom 



that the sum of (—1)^ for the configurations on 5*2 van- 
ishes as soon as at least one site on iS'2 can be both empty 
and occupied. Consequently, the only non-zero contribu- 
tion comes from the configurations such that at least one 
of the adjacent sites on 5*1 is occupied. There are only two 
such configurations: 

|a) = ••••□ o •□ o •□ o •□ o •□ o •□ o •□ o . . . 

I7) = • • • o □ • oD • on • oD • oD • oD • on • . . . (3) 

where the square represents an empty site on 6*2 . The final 
step is to sum (— f)''^ for these two configurations, and 
since both configurations have f = L/3 = j, we find that 
the Witten index is W = 2{—iy. Note that this agrees 
with the result we obtained for the 6-site chain. 

To find the exact number of ground states we compute 
the cohomology by using a spectral sequence. A useful the- 
orem is the 'tic-tac-toe' lemma of [8,. This says that under 
certain conditions, the cohomology Hq for Q ^ Q1 + Q2 is 
the same as the cohomology of Qi acting on the cohomol- 
ogy of Q2- In an equation, Hq — Hq-^ {HQ2) = H12, where 
Qi and Q2 act on different sublattices and 82- We find 
H12 by first fixing the configuration on all sites of the sub- 
lattice Si , and computing the cohomology Hq^ . Then one 
computes the cohomology of Qij acting not on the full 
space of states, but only on the classes in Hq^. A suffi- 
cient condition for the lemma to hold is that all non-trivial 
elements of H12 have the same /2 (the fermion-number on 
5*2). For the periodic chain with L = 3j we choose the 
sublattice as before. Now consider a single site on 82- If 
both of the adjacent sites are empty, Hq^ is trivial: Q2 
acting on the empty site does not vanish, while the filled 
site is Q2 acting on the empty site. Thus Hq^ is non-trivial 
only when every site on 5*2 is forced to be empty by be- 
ing adjacent to an occupied site. The elements of Hq^ are 
just the two states |a) and I7) pictured above in eq. ([3]). 
Both states \a) and I7) belong to Hi2'- they are closed be- 
cause Qi\a) = Q1I7) = 0, and not exact because there are 
no elements of Hq^ with /i = / — 1. By the tic-tac-toe 
lemma, there must be precisely two different cohomology 
classes in Hq, and therefore exactly two ground states 
with / = L/3. Applying the same arguments to the peri- 
odic chain with 3/ ± 1 sites and to the open chain yields in 
all cases exactly one E = Q ground state, except in open 
chains with 3/ -I- 1 sites, where there are none [5]. 

The supersymmctric model on the chain can be solved 
exactly through a Bethe Ansatz . In the continuum limit 
one can derive the thermodynamic Bethe ansatz equa- 
tions. The model has the same thermodynamic equations 
as the XXZ chain at Z\ = —1/2, so the two models coincide 
(the mapping can be found in [5]). The continuum limit 
of the XXZ chain is described by the massless Thirring 
model [9], or equivalently a free massless boson with 
action fTU] 

8=ljdxdt [{dtl^f ~ [d^l^f] . 

The continuum limit of the A = —1/2 model has g = 2/3; 
this is the simplest field theory with N = (2, 2) supercon- 
formal symmetry fTO^. The (2, 2) means that there are two 



left and two right-moving supersymmetries: in the contin- 
uum limit the fermion decomposes into left- and right- 
moving components over the Fermi sea. So the system is 
critical and in the continuum limit it is described by a 
superconformal field theory. 

A final note we would like to make here, is that the 
supersymmctric model on the chain recently emerged in 
a special limit of a large- A'^ supersymmctric matrix model 
|llj . This is yet another interesting connection, worthy of 
further investigation. 



3 Beyond ID: heuristics and methodology 

We have seen that in the one dimensional case the Hamil- 
tonian favors a configuration where the hard-core fermions 
sit three sites apart on average. This heuristic picture can 
be extended beyond ID. For convenience we restate the 
Hamiltonian in its general form 

i j next to i i 

The second term in the Hamiltonian gives a positive con- 
tribution to the energy for every site that has all neigh- 
boring sites empty, regardless of whether the site itself 
is occupied or empty. For every hard-core fermion it will 
give a contribution of +1, since by definition it has its 
neighboring sites empty. So the contribution of this term 
is minimized by blocking as many sites as possible with as 
few fermions as possible. Roughly speaking, this criterion 
leads to configurations where all fermions are three sites 
apart. We call this the '3-rule'. In the following sections 
we shall see this 3-rule in action, and demonstrate how it 
leads to superfrustration. 

In the previous section we have developed some tools 
to study our model on different lattices. The Witten index 
gives us a lower bound on the total number of supersym- 
mctric ground states. In some cases one can find a recur- 
sion relation or even a closed form for the Witten index 
as a function of the system size. The growth behavior of 
this function gives a lower bound to the growth behavior 
of the ground state entropy. 

If the Witten index grows exponentially with the sys- 
tem size, we have an extensive ground state entropy. In 
Sect. 14.11 we present an example where this entropy is 
known in closed form. Numerical studies of the Witten in- 
dex on two-dimensional lattices [12] have revealed that for 
generic lattices the Witten index indeed grows exponen- 
tially with the system size (see Sect. 14.3] for an example). 
The square lattice is an exception, since there the Wit- 
ten index grows exponentially with the perimeter of the 
system. We shall touch upon some features of the square 
lattice in Sect. 14.41 

Further insight in the ground state structure can be ob- 
tained from the cohomology method. This gives the exact 
number of ground states and for each of them the number 
of fermions. A remarkable feature of our model that has 
been revealed by cohomology studies is that ground states 
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typically occur at different fermion- numbers, or equiva- 
lently at different filling fractions. This implies that the 
Witten index will typically underestimate the actual num- 
ber of ground states. More remarkably, however, this also 
implies that one can add a particle to the system or ex- 
tract a particle from the system within a certain window 
of filling fractions without paying any energy. 

Finally, it has in many cases turned out to be pos- 
sible to characterize supersymmetric ground states with 
the help of an 'effective geometric picture'. In this, one es- 
tablishes a (almost) 1-1 correspondence between quantum 
ground states and geometric configurations such as cover- 
ings of the lattice by dimers or tiles of specific dimensions. 
Examples are dimer coverings for the case of the martini 
lattice (Sect. I4.ip and rhombus tilings for the 2D square 
lattice (Sect. 133) • The geometric picture is related to the 
heuristic 3-rule, but much more robust. It has been pio- 
neered in [3| and in a remarkable series of mathematical 
papers by J. Jonsson pTJ . 

The fact that the ground states of strongly correlated 
quantum fermion models can be characterized by geomet- 
ric means is quite deep and at this time not fully under- 
stood. It suggests that further properties of these models 
(such as the excited state spectra) are tractable by similar 
means, which opens most interesting perspectives. 

It the next section we present a gallery of examples 
of 2D lattices, and indicate what has been revealed about 
their ground state structure using the various approaches 
described in this section. A more systematic account is 
forthcoming [T3]. 

4 Beyond ID: examples 
4.1 Martini lattice 

The martini lattice (see Fig. [J) is an example of a two 
dimensional lattice, where the cohomology can be com- 
puted relatively easily [3] . The method is strongly related 
to the one used to compute the cohomology for the chain. 
The computation proves that the ground state entropy is 
an extensive quantity and we find a closed expression for 
the ground state entropy per site. We shall see that the 
3-rule is not violated in this case. This is related to the 
fact that the martini lattice, due to its structure, nicely 
accommodates the 3-rule. Lattices with a higher coordi- 
nation number usually do not have this property and con- 
sequently allow for a window of filling fraction for the 
supersymmetric ground states. 

The martini lattice is formed by replacing every other 
site on a hexagonal lattice with a triangle. To find the 
ground states, take Si to be the sites on the triangles, 
and ^2 to be the remaining sites. As with the chain, Hq^ 
vanishes unless every site in S2 is adjacent to an occupied 
site on some triangle. The non-trivial elements of Hq^ 
therefore must have precisely one particle per triangle, 
each adjacent to a different site on 82- This is because a 
triangle can have at most one particle on it, and (with ap- 
propriate boundary conditions) there are the same number 
of triangles as there are sites on S'2. A typical element of 




Fig. 2. Hard-core fermions on the martini lattice (left) and on 
the kagome ladder (right). 

Hq^ is shown in Fig. [21 One can think of these as 'dimer' 
configurations on the original honeycomb lattice, where 
the dimer stretches from the site replaced by the triangle 
to the adjacent non-triangle site. Each close-packed hard- 
core dimer configuration is in H12, and by the tic-tac-toe 
lemma, it corresponds to a ground state. The number of 
such ground states e^'-^^ is therefore equal to the number 
of such dimer coverings of the honeycomb lattice, which 
for large L is [15] 

S'rs 1 

— = - de ln[2 cos e] = 0.16153 . . . 

The frustration here clearly arises because there are many 
ways of satisfying the 3-rule. 

4.2 Kagome ladder 

In this section we consider the kagome ladder (see Fig. [2]) 
as an illustration of a case where we can compute the co- 
homology exactly, but where the 3-rule is not very helpful. 
We find a closed expression for the partition function and 
a window of filling fraction for the supersymmetric ground 
states, which can both be interpreted in terms of tilings. 

Computing the cohomology in this case is a bit more 
involved due to two things: First, in the previous examples 
we could always choose the sublattice S'2 such that it con- 
sisted of disconnected sites, which by themselves have zero 
cohomology. For the kagome ladder the convenient choice 
for the sublattice S'2 is less trivial. The second complica- 
tion arises because not all elements of the cohomology H12 
will have the same fermion-number, which was a sufficient 
condition for the tic-tac-toe lemma to hold. 

We compute the cohomology step by step. For each 
step there is a supporting picture in Fig. [3| Step 1 is 
to map the kagome ladder with hard-core fermions to a 
square ladder with hard-core dimers. This mapping is one- 
on-one. In Step 2 wc define the sublattices Si and S2. 
Step 3 is to note that the cohomology of one square plus 
one additional edge vanishes. To do so, first note that if 
site 3 and 4 are both empty, we get zero cohomology due 
to site 5 which can now be both empty and occupied. The 
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Step 1: 



Step 2: 



Step 3: 



Step 4: 
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Fig. 3. Step-by-step computation of the cohomology of the 
kagome ladder. 



remaining four configurations are easily shown to be ei- 
ther Q of something (exact) or not in the kernel of Q (not 
closed) . Step 4 is to build up the ladder by consecutively 
adding blocks with 3 rungs (9 edges) to the ladder. From 
step 3 we now conclude that there are only two allowed 
configurations for the additional ^i-sites: they must be ei- 
ther both empty or both occupied, since if just one of them 
is occupied the remaining configurations on the additional 
S'2-sites are exactly the ones of step 3. A simple computa- 
tion shows that if both sites on Si are occupied, Hq^ has 
one non-trivial element with one dimer and if both sites 
on 5*1 are empty, Hq^ has two non-trivial elements, both 
with two dimers. 

Now it is important to note that on the three addi- 
tional rungs we find three non-trivial elements of the H12, 
but two of them have /2 = 2 and one has /2 = 1. It 
can be shown that all three indeed belong to Hq by go- 
ing through the tic-tac-toe lemma step by step. This is a 
tedious computation, but it can be done. 

Finally, we find the cohomology of the kagome ladder 
with open boundary conditions of length n, which corre- 
sponds to a ladder with n + 1 rungs and 3n -I- 1 edges in 
total, by recursively adding rungs to the system. We thus 
obtain a recursion relation for the ground-state generat- 
ing function P„(z) = trGs(-z^), which gives the Witten 
index for z = — 1 and the total number of ground states 
for z = 1: 

Pn+3{z) = 2z2p„(z) + Z^Pn-l{z), 

with Pq = 0, Pi = z, P2 = 2z'^, P3 = z^. Instead of draw- 
ing conclusion from here, let us picture the above in terms 
of tiles. From step 4 we conclude that we can cover the 
ladder with three tiles, two of size 9 (i.e. 9 edges) contain- 
ing 2 dimers and one of size 12 containing 3 dimers. From 
this picture we obtain the same recursion relation pro- 
vided that we allow four initial tiles corresponding to the 
initial conditions of the recursion relation above. Further- 




Fig. 4. The M X N triangular lattice has periodic boundary 
conditions along the directions of the two arrows. 

Table 1. Witten Index for the M x N triangular lattice. 
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5237 


11 


1 


67 


1 


859 


1651 


12607 


32418 


12 


1 


-47 


130 


-1295 


-589 


-26006 


-152697 


13 


1 


-181 


1 


-77 


-1949 


67523 


330331 


14 


1 


-87 


-257 


3641 


12611 


-139935 


-235717 


15 


1 


275 


-2 


-8053 


-32664 


272486 


-1184714 



more, we can see directly that the window of filling frac- 
tion of the tiles runs from 2/9 to 1/4. Using the recursion 
relation, we find that the ground-state entropy is set by 
the largest solution Amax of the characteristic polynomial 
A4 _ 2A - 1 = 0, giving Sgs/L = (In A,„ax)/3 = 0.1110 ... . 



4.3 2D triangular lattice 

The ground state structure of the supersymmetric model 
on the 2D triangular lattice is not fully understood. Never- 
theless, it is clear that ground states occur in a finite win- 
dow of filling fractions v = f / L and that there is extensive 
ground state entropy. These features seem to be generic 
for 2D lattices, as they have been observed for many ex- 
amples such as hexagonal, kagome, etc. (2D square being 
an important exception) |12j . 

In Tab. [T] we show the Witten indices for the M x 
N triangular lattice, with periodic boundary conditions 
applied along two axes of the lattice (see Fig. [4|). The 
exponential growth of the index is clear from the table. 
To quantify the growth behavior, one may determine the 
largest eigenvalue Atv of the row-to-row transfer matrix 
for the Witten index on size M x N. This gives 



A 



N 



\WmM ~ (Ajv)" + (Xn) 
|A|-1.14, arg(A) - 0.18(7r) 
leading to a ground state entropy per site of 
Sgs 



A 



N 



> 



1 



MN - MN 



loglM^M, 



N 



log|A| - 0.13 . 



(4) 



(5) 
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Table 2. Witten Index for M x A'' square lattice. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


1 


-1 


1 


3 


1 


-1 


1 


3 


1 


-1 


1 


3 


3 


1 


1 


4 


1 


1 


4 


1 


1 


4 


1 


1 


4 


4 


1 


3 


1 


7 


1 


3 


1 


7 


1 


3 


1 


7 


5 


1 


1 


1 


1 


-9 


1 


1 


1 


1 


11 


1 


1 


6 


1 


-1 


4 


3 


1 


14 


1 


3 


4 


-1 


1 


18 


7 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


8 


1 


3 


1 


7 


1 


3 


1 


7 


1 


43 


1 


7 


9 


1 


1 


4 


1 


1 


4 


1 


1 


40 


1 


1 


4 


10 


1 


-1 


1 


3 


11 


-1 


1 


43 


1 


9 


1 


3 


11 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


12 


1 


3 


4 


7 


1 


18 


1 


7 


4 


3 


1 


166 


13 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


14 


1 


-1 


1 


3 


1 


-1 


-27 


3 


1 


69 


1 


3 


15 


1 


1 


4 


1 


-9 


4 


1 


1 


4 


11 


1 


4 



The argument of A indicates that the asymptotic behavior 
of the index is dominated by configurations with fiUing 
fraction around v = 0.18. 

In a most interesting mathematical analysis 13J, Jon- 
sson has shown that for a sufficiently large triangular lat- 
tice (with open BC) ground states occur for all rational 
numbers in the range 



1/7 <v< 1/5 



(6) 



His analysis is based on an effective geometric picture in- 
volving so-called cross-cycles, which however is less ex- 
plicit than the one that has been worked out for the case 
of the 2D square lattice (see below) . There is a clear chal- 
lenge to develop this picture further to the point that the 
growth behavior of the Witten index, and of the number 
of ground states, can be given in closed form. 



4.4 2D square lattice 

Numerical studies [16] of the Witten index of the square 
lattice revealed a very different behavior (see Tab. [2]). At 
first glance one notices that it does not grow exponentially 
with the system size. In fact, more detailed investigation 
of these studies led to two conjectures [1^ for which a 
proof was found by Jonsson [T^. We state one of these 
results here: 

for an M x N square lattice with periodic boundary 
conditions in both directions, W — 1 when M and 
TV are coprime. 

Extending this work, Jonsson found a general expres- 
sion for the Witten index Wu,v of hard-core fermions on 
the square lattice with periodic boundary conditions given 
by the vectors u = (ui,U2) and v = (wi,W2)- The M x N 
square lattice is now a specific case with u = [M, 0) and 
V = (0, N) (for an extension of this work to other fami- 
lies of grid graphs see [IT])- A crucial step in [13] is the 
introduction of rhombus tilings of the square lattice. It is 
shown that the trace in the Witten index can be restricted 
to configurations that can be mapped to coverings of the 
plane with the four rhombi or tiles shown in Fig. (5] Note 
that the sides of these rhombi, which connect the hard- 
core fermions, are in agreement with the heuristic 3-rule. 
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Fig. 5. Tilings of the 2D square lattice. Above: the four dif- 
ferent rhombi. Below: mapping between tiles and hard-core 
fermions. 



Furthermore, two of the rhombi have area 4, whereas the 
other two have area 5. A covering with either of these 
rhombi alone thus corresponds to a filling fraction of 1 /4 
or 1/5, respectively. 

To state Jonsson's results for the Witten index we in- 
troduce the following notations. We denote by Ru,v the 
family of tilings of the plane with boundary conditions 
given hy u = (mi, U2) and v = (wi, ^2). Furthermore 
and are the number of tilings of this plane with an 

even and an odd number of tiles, respectively. Finally, we 
define 

2 if d = 3fc, with k integer 
— 1 otherwise. 



(7) 



The expression for the Witten index then reads [13] 

Wu,v = -{-l)%0,, + \R+J - \R-J, (8) 

where d = gcd(Mi — U2, vi — V2) and d* = gcd(wi -f U2, vi + 
W2). It can be shown that the Witten index grows exponen- 
tially with the linear size (not the area) of the 2D lattice. 
Detailed results for the case of diagonal boundary condi- 
tions have been given in [18] . Further studies of the Witten 
index transfermatrix for the square lattice with diagonal 
and free boundary conditions by R.J. Baxter [19] have led 
to an additional set of conjectures. 

As we already mentioned in Sect. |3l the geometric 
picture in terms of tilings is useful beyond the compu- 
tation of the Witten index. It also provides a way to de- 
termine a window of filling fractions in which supersym- 
metric ground states can be found. The result is [13] that 
for large enough square lattices (with open BC) ground 
states exist for all rational fillings in the range 



1/5 < < 1/4 



(9) 



While it is clear that the effective geometric picture 
in terms of rhombus tilings goes a long way characteriz- 
ing the supcrsymmetric ground states, it has until now 
failed to give complete results for the ground state parti- 
tion sum for the supcrsymmetric model on the 2D square 
lattice. For this issue, and for many others, it is impor- 
tant that the results presented in this section hold for the 



8 



Liza Huijse, Kareljan Schoutens: Superfrustration of charge degrees of freedom 



square lattice with any kind of periodic boundary con- 
ditions. This also includes seini-2D lattices, i.e. various 
ladders and even the ID chain. These lattices are a good 
arena to further investigate properties of the supersym- 
metric fermion models, both analytically and numerically 



5 Conclusion 

The analysis of strongly correlated fermions on lattices in 
dimension D > 1 is a notoriously difficult problem, for 
which very few exact results have been obtained. At the 
same time, the problem is highly relevant, as it holds the 
key to the behavior of correlated electrons in quasi-2D ma- 
terials. We have here presented various exact results for 
the ground state structure of a fermion lattice model with 
an exact supersymmetry. In particular, we have demon- 
strated the remarkable feature of superfrustration, which 
this model possesses on generic 2D lattices. 

In our discussion of the various examples of superfrus- 
tration, we mostly focused on specifying the number of su- 
persymmetric ground states and the fermion number (or 
filling fraction) where they occur. Clearly, one would like 
to understand better various properties of these states, as 
well as the quantum phases they give rise to when param- 
eters are perturbed away from the supersymmetric point. 

The supersymmetric ground states on the ID chain are 
quantum critical and as such described by a superconfor- 
mal field theory. For a more general class of supersymmet- 
ric ID models [5] (where the nearest neighbor exclusion 
rule is softened) the situation is akin to that of higher- 
S spin chains: the models are gapped but go critical if 
interaction parameters are tuned to specific values. For 
the 2D models presented here, the issue of quantum crit- 
icality is under investigation [14]. While supersymmetry 
alone certainly does not imply quantum criticality, it is 
clear that the balancing between kinetic and interaction 
terms that is implied by supersymmetry steers one into 
regions of parameter space where charge order and Fermi 
liquid behavior compete. 
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